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Abstract 

We analyze the signatures of inflationary models that are coupled to strongly interacting field the- 
ories, a basic class of multifield models also motivated by their role in providing dynamically small 
scales. Near the squeezed limit of the bispectrum, we find a simple scaling behavior determined 
by operator dimensions, which are constrained by the appropriate unitarity bounds. Specifically, 
we analyze two simple and calculable classes of examples: conformal field theories (CFTs), and 
largc-N CFTs deformed by relevant time-dependent double-trace operators. Together these two 
classes of examples exhibit a wide range of scalings and shapes of the bispectrum, including 
nearly equilateral, orthogonal and local non-Gaussianity in different regimes. Along the way, we 
compare and contrast the shape and amplitude with previous results on weakly coupled fields 
coupled to inflation. This signature provides a precision test for strongly coupled sectors cou- 
pled to inflation via irrelevant operators suppressed by a high mass scale up to ~ 10^ times the 
inflationary Hubble scale. 
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1 Introduction and summary 



In inflationary cosmology, the quantum fluctuations of the inflaton, and possibly other fields, get 
imprinted on the power spectrum and higher point correlations accessible in the CMB and large- 
scale structure. This basic idea goes back to the origins of the subject, and has been explored 
in many illustrative examples. Recently there has been progress toward a more systematic un- 
derstanding of the observables and their implications, with input both from bottom up effective 
field theory and from top down UV complete mechanisms for inflation. 

In this paper, we mix in a basic class of strongly-coupled fields - field theories which are 
conformal at the Hubble scale and their time-dependent deformations - and compute their con- 
tribution to the perturbation spectrum and non-Gaussianity. We find characteristic scaling be- 
havior depending on the anomalous dimensions of operators, and compute in detail the shape 
and amplitude of the non-Gaussianity. The results depend on some basic properties of quantum 
field theory, such as unitarity bounds and crossing relations. 

Aside from theoretical interest, this analysis is motivated by several considerations. First, 
it is a canonical and calculable regime of field theory which belongs in a systematic analysis 
of multifield signatures and their significance, particularly given the observational handles on 
non-Gaussianity. Relatedly, as we will see, these strongly coupled sectors exhibit an interesting 
partial degeneracy with weakly interacting fields in their predictions for the simplest practical 
observables. Finally, from the top down, strongly coupled sectors play a useful role in model- 
building, producing naturally small scales in the effective action for the inflaton, and fields from 
this sector may participate in the perturbations. More generally, we can use forthcoming non- 
Gaussianity data (whether a constraint or detection) to provide a precision test of additional field 
theory sectors which couple to inflation via higher dimension operators. 

To put this in context, one basic question probed by observations is the number of fields 
which participate in generating the correlation functions of the curvature perturbations. This is 
a different question from the question of how many fields are involved in the underlying mechanism 
producing the inflationary background. One can separate multi-field theories into two classes, 
depending on whether the additional fields themselves acquire scale invariant perturbations that 
affect the curvature or the isocurvature fluctuations. In the case we will be interested in here, the 
additional flelds instead affect the density perturbations through their couplings to the inflaton. 
Though not directly observable, these additional fields cannot usefully be integrated out, as doing 
so would lead to a non-local Lagrangian. Progress toward a general treatment of both types of 
theories can be found in [1][2], where the effects of additional sectors are packaged in terms of 
their correlation functions. 

A useful general theorem known as the consistency condition establishes that a single-field 
theory of perturbations cannot generate non-Gaussianity in the 'squeezed' configuration, with one 
mode much longer than the others [3j. Conversely, it has been long known that non-Gaussianity 
peaked in the squeezed configuration can arise in the presence of additional light fields [4J. For 
massless weakly interacting fields coupled to the inflaton, the three point function of scalar 
perturbations behaves as 



(CfciCfcaCfca) 



fmA^ 



(27r)3,5(=')(ki+k2 + k3) 




2 



where ~ 10 ^ is the ampUtude of the scalar perturbation. Moreover, interesting power-law 
deviations from this shape of non-Gaussianity arise for fields of nonzero mass m [5] 

iCkMks) = B{h,k2,h) {2Trf6^^Hk,+-k2 + 'k3) (1.2) 

^^1) (2.)V')(k,+k, + k,) (1,3) 

This scenario, known as quasi single-field inflation, is a feature of weakly coupled theories with 
supersymmetry broken at the Hubble scale [6j. 

In the present work, we will analyze two simple cases where the additional fields affecting the 
density perturbations are strongly interacting, using conformal symmetry and related methods 
to control the calculations. For our first class of examples, we will study inflation coupled to a 
conformally coupled CFT, for which the conformal dimension A of an operator O plays the role 



of the mass-dependent exponent in (1.3). For our second class of examples we will consider a 
simple time-dependent flow away from a large-N CFT, which gives unitary theories realizing a 
larger range of exponents, including local /nl. 

For our first class of examples, a conformally coupled CFT linearly mixes with inflaton per- 
turbations , giving rise to a bispectrum of the form 

^^^^^^ ^xf^V (A>2). 

For this case, the standard unitarity bound implies A > 1. As a result, the shape is peaked 
at equilateral/flattened triangles in momentum space, with the scaling behavior (1.4) as one 
approaches the squeezed limit determined by the dimension of the most relevant operator that 
couples in. Moreoever, we will find a range of exponents arising from dimensions 3/2 < A < 2 
which cannot be obtained in quasi-single field models, producing a scale-dependent bias which 
could distinguish them. For the range of exponents where the two are degenerate, it is intriguing 
that our intrinsically gapless CFT fields behave like massive weakly coupled fields with respect 
to the squeezed limit; the origin of this effect is the redshifting of conformally rescaled correlators 
of our conformally coupled operators. The conformal coupling to curvature which goes into this 
analysis is a special choice, and we expect a wider range of behaviors in the presence of more 
general curvature couplings. 



Indeed, our second class of examples will give rise to the same scaling as (1.4) while allowing 



for A < 1. These examples make use of the fact that time-dependent couplings in quantum field 
theory can strongly affect infrared physics and shift unitarity bounds, as studied recently in [^. 
This is important in the present context since time dependent couplings can arise very easily 
via couplings of the rolling infiaton field to other sectors such as a CFT. In the examples in [,7], 
time dependent couplings can introduce fiows between a unitary CFT containing an operator O 



^Note that our CFT lives in the four-dimensional spacetime, and is not to be confused with conjectural lower- 
dimensional holographic duals for de Sitter, 
^up to logarithmic factors. 



3 



of dimension and an infrared theory with two-point correlators fahing off like 1 / distance^^- 
times powers of the time-dependent couphng, where A_ = 4 — A-|_. In particular a theory with 
a marginal scalar operator with A-i- ~ 4 can flow in this way to a theory with A_ ^ 1, giving 
nearly local non-Gaussianity, not suppressed in the squeezed limit by any additional powers of 

fcl/fc2. 

Finally, let us discuss the amplitude Jml of the bispectrum. As we will see, this can be 
substantial, and will give us sensitivity to higher dimension couplings of the inflationary sector 
to other fields. To give a rough illustration of this last point, consider, for example, a dimension 
A = 2 operator O2 in a CFT, coupled to the inflaton (j) via the dimension six operator 

.'.^<^ . (..5) 

Let us evaluate one factor of dtj) on the backround rolling scalar field, which in slow-roll inflation 
given in terms of the inflationary Hubble scale by 0o ~ H"^ / ^C, ~ lO^i?^. This gives us a 
linear mixing J A^^{H/M^)'^6(j)0 between the canonical perturbation 6(f) and oj^ The CFT 
three-point function combined with three of these mixing interactions generates a contribution 
to the three-point function of the inflationary perturbations. This leads to an amplitude f^L 
behaving parametrically like 



/^L~CA^-V(2)(^) (1.6) 



where C is the amplitude of the CFT three-point function (2.1), and the function /(A), plotted 
below in Fig. 6, gives a substantial numerical factor]^ Since current measurements are projected 
to be sensitive to Jnl ^ 10, they provide a precision test of this higher dimension coupling 
up to a value of M^: ~ lO'^H. This high mass scale (relative to Hubble) is generally below the 
Planck scale (Mpi), but can exceed the GUT scale in simple examples such as chaotic inflation. 
In string-theoretic ultraviolet completions of inflation, higher dimension operators suppressed by 
a scale Af=K <C Mpi of order the Kaluza-Klein or string scale may be probedj^ That is, while an 
observation of non-Gaussianity may be explained by contributions from additional flelds (weakly 
or strongly coupled), a null result would conversely provide a precision constraint on very high- 
energy physics. Of course sensitivity to high energy physics arises already in single field inflation 
(see e.g. [12 ^ 1131 [T^ I15| ). A feature of the present case is that observations will constrain hidden 



sectors of additional fields coupled to infiation via higher dimension operators such as (1.5). 
Similar constraints arise in other models involving a linear mixing with the infiaton (e.g. [5l I16j). 

Both of our main examples are just calculable examples of a wider point: additional fields 
active during infiation may include strongly coupled sectors with characteristic signatures. It 

''For the purposes of finding a simple estimate of our sensitivity to high dimension couplings, we will tune away 



the relevant perturbation of the CFT that we get from ( 1.5 1 if we evaluate both factors of d(j} on the rolling inflaton 



background. In the detailed examples to be discussed in the main body of this paper and appendix A, we will find 
similar results without such tuning. 

''Below we will describe our normalization conventions which go into this. 

^In several of the UV completions of inflation explored extensively in string theory, reviewed for example in [S], 
strongly coupled sectors play a useful role in producing dynamically small scales [9] . Operators suppressed by the 
Kaluza-Klein scale are ubiquitous there [10]. (See [11] for some more recent examples in which interacting field 
theory plays a role in the inflationary mechanism.) 
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will be interesting to analyze this more generally and systematically, and also to incorporate the 
particular couplings arising in complete models of inflation which involve couplings to strong 
dynamics to obtain their specific multifield signatures. 

2 Conformally Coupled Examples 

In this section, we will consider strongly coupled theories which behave like a CFT near the 
inflationary Hubble scale. In flat spacetime, CFTs are one of the best studied classes of interacting 
field theories. Due to the high degree of symmetry, much is known about the spectrum of 
operators and their correlations functions. Because de Sitter space is conformally flat, a CFT 
can be coupled to gravity such that correlation functions in de Sitter space preserve the flat 
space results up to an overall rescaling. This choice is related to the choice of curvature couplings 
involving operators in the CFT (or equivalently to the choice of improvement terms for the stress 
tensor that we couple to gravity). In this section, we will assume these couplings are chosen to 
preserve conformal invar iance. 

2.1 Setup 

We are interested in the possibility that a CFT is weakly coupled to the inflaton and its perturba- 
tion. As a result, the CFT will influence the correlation functions of the curvature perturbation 
we observe at late times. For this purpose, it is not necessary to specify the underlying dynamics 
leading to inflation; instead it is most convenient to work directly in terms of the perturbations 
7r(t, x) as in the effective field theory treatment developed in |14j. Indeed, the effects we wil 
compute could accompany a wide variety of underlying inflationary mechanisms. 

In the absence of inflation, a CFT is described by a list of local, primary operators 0^'-'^'^' (x, t) 
and their correlation functions, where € (|, |) is the spin and Aj is the dimension of the 

i-th operator. For simplicity we will focus on a single scalar operator (D{x,t) with dimension A. 
In flat space, the two- and three-point functions of 0(x, t) are fixed up to a constant. Specifically, 
the correlation functions in Euclidean signature take the form 

(O(x,t)O(0)) = ^2^^, (2.1) 
(0(xi,ti)0(x2,t2)0(x3,t3)) = A ^^-^ A' (2-2) 

1^12 + ^12! ^ 1^23 + ^23! ^ 1^31 + *3ll ^ 

where is a constant. The normalization of the two-point function 

is a convention [3 

The metric of de Sitter space is given by 

ds^ = a{T)^{-dT^ + dx^) . (2.3) 

^In the limit A = 1, the CFT becomes free. In this case, the convention for the two-point function differs from 
the one of a free scalar field a factor of 47r^. This explains why our dimensionless functions have sometimes large 
numerical values. 
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where a(T) = —1/{Ht). This metric is conformaUy flat in terms of the conformal time r. 
Therefore, we get the CFT correlation functions in de Sitter space from the flat space result by 
replacing t — )• r and a Weyl transformation of 0{x,t) — )• a(r)~^0(x, r). In particular, the two 
and three point functions are given by (in Euclidean time): 

(O(r.x)O(/.x0> = , (2.4) 

[(r - r')^ + {x- x')^] 

(0(xi,ri)0(x2,r2)0(x3,T3)) = , 9|A/2|^2 , ^2|A/2|^2 , ^2 |A/2 (^.S) 

1-^12 "'"^21 1-^23 """^231 1-^31 "T^SlI 

In particular, these correlation functions redshift at late times much like those of a weakly inter- 
acting massive field. These strongly coupled correlators are different from those that arise in the 
weakly coupled case [5] [6] ; the position space two point function for massive fields is a nontrivial 



hypergeometric function, distinct from the simpler function (2.4) except for the special case of a 
conformally coupled free scalar. 

We will now couple the inflationary perturbation vr to this CFT via the interaction Hamilto- 
nian 

^int = \f?-^M^^\H\^'^{2ii - d^,7Td>'TT)0 + \M^^\H\fl-^ {-2ir + d^-Kd'^iifO . (2.6) 

Here and elsewhere, / = —Ht^ is a derivative with respect to FRW time t = — \og{—HT)/H. 
The scalar vr ~ 5t is not canonically normalized; it is related to the canonically normalizecj^ 
perturbation tTc via vr = tTc/ \J 2M^^H. It is also related to the conventionally normalized scalar 
perturbation C via C = —Hn (at linear order). We also note here that at freezeout, the corre- 
sponding field amplitudes which will enter into the calculations below are ( ~ ~ 10~^ and 
TTc ^ H. In the whole paper we will neglect the mixing with gravity, as it will give subleading 
corrections. 



In a slow-roll model like (1.5), the parameters fi and jl may be related to (p and some higher 



scale, (e.g. /j, = (f)/M^, with ^ (j) when A = 1). However, there may be many other 



UV completions that also give rise to (2.6) where these scales have different origins. For this 
reason, we will work directly with fi and jl throughout, as they are the parameters relevant to 
the phenomenology. 

When A < 2 the leading contribution from the first term is a relevant deformation and 
therefore is perturbative when fj, <^ H. On the other hand, the second term is irrelevant for 
all dimensions consistent with unitarity and is therefore perturbative the jl ^ H. A priori, 
the second term may or may not contribute significantly to the bispectrum. We will therefore 
consider the two cases separately in section [2^ and [2?5l For the special case of A = 2, we should 
replace (/x/i^)^^^ with a dimensionless coupling A: A = limA^>2(/w/-f^)^~^i with log{fi/H) — )• 
(A — 1)/(A — 2). This replacement should be unambiguous so we will not do it explicitly. 

In Appendix A, we analyze the radiative stability of this setup. One result of that analysis is 
that under appropriate conditions the term ~ J m^'^O generated by vr loops satisfies m H, 
meaning that even for relevant operators (A < 4) we do not generate a flow away from the CFT 
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Here we are assuming unit speed of sound, Cs ~ 1. For general sound speed, the canonically normalized field 



is given by tt = CsiVc/ \/'2AKH 
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over the scales of interest. This analysis is self-contained up to a scale A which can be ^ H, 
leading to precision tests of higher dimension operators as anticipated in the introduction. 

2.2 Calculating in-in Correlators in Euclidean Signature 

Throughout the paper, we will be interested in calculating in-in correlations functions of C 
evaluated at equal times. These can be computed perturbatively, using the interaction picture 
fields, stating from [T7j, 

{fexp[i [ i7int(T)a(T)dr]('Cint(ki,ro)..Cint(kn,ro))rexp[-i j i?int(T)a(r)dr]) , 

J-oo{l+it) ^ ' J-oo{l-it) 

(2.7) 

where we have assumed the Bunch-Davies vacuum and i7int(T) = / d^x a'^{T)'HintiT, x). In order 
to simplify these calculations, we will follow the strategy suggested in [TH]. The basic idea is that, 
after rotating the conformal time integrals via r — )• ztir^; + tq, the in-in correlation becomes a 
(Euclidean) anti-time-ordered correlation function. 




Figure 1: The analytic continuation of the contour in conformal time (r) from Lorentzian signature (blue) to 
Euclidean signature (red). We calculate the correlation functions for operators at tq < 0. Our expressions involve 
branch cuts only when Rer > 0, which ensures this continuation is well defined. 



Starting from the ie prescription in (2.7), which projects onto the interacting vacuum, the 
time and anti-time ordered exponentials should be rotated to opposite values of the Euclidean 
time. After doing this, we obtain 

{Texp[i[" ifint(T)a(r)(ir](Cint(ki,To)..Cint(k„,ro))rexp[-i / ° Fint(r)a(r)dT]) . (2.8) 
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This we recognize as simply the anti-time ordered correlation function in Euclidean time 



(f (Cint(ki,ro)..Cint(kn,ro)exp[f / /7i„t (r)a(T)dr] J) (2.9) 

= (rKint(ki,ro)..Cint(k„,ro)exp[- / iJint(iT£; + ro)a(fr£' + ro)(iT£;] j) . (2.10) 

J — CO 

Provided we use the anti-time ordered, Euclidean Green's function which we will compute mo- 
mentarily, the operator ordering is automatic. This analytic continuation is illustrated in figure 

m 

The anti-time-ordered Green's functions for any two operators can be written as 

(tOi(r)02(r')> ^ e{T'E - te){0i{t)02{t')) + 0{te - t'e){02{t')0i{t)) , (2.11) 

where 9{x) is the Heaviside step function. Following the usual quantization of a scalar in de 

+ ^fe('^)^k with [ok, a^, 



Sitter space, we write the operator vTk = '^k{T)0't + ^fe(''")ok with [a^, at,] = (27r)^(5(k + k') and 



^ ^ (1 - ikr) H (1 + ikr) 

^'^^^ = 2M,my^ k^l- ' ' "^^"^ = 2Mp,|if|V2 ^3/2 - • (2-12) 

(As usual, we conjugate tt to obtain vr* in our original Lorentzian signature calculation, before 
deforming our contour to lie along r^;.) Notice that external factors of 7rk(To) only have nonzero 
contractions with factors of tTc in T^mti so that we will only need 

(r(^e(irE+To,ki)7re(ro,k2))) = -^(l+iA:iro^)(-iri5-ro)2e-'=il-^l(27r)35(ki+k2) . (2.13) 
V / lk\ \te\ 

The appearance of the absolute values in the above expression does not obstruct the analytic 



continuation we performed in writing (2.8); the only non-analytic behavior appears at the location 



of the operator insertions, r = tq, which is fixed in our continuation. Alternatively, one can choose 



the domains of integration in (2.7) to be manifestly time-ordered, which ensures the integrands 
are analytic in r. 

As we discussed in the previous section, the GET correlation functions in Euclidean signature 
may be be taken to be anti-time ordered. Therefore, when evaluating correlation functions we 



can use equations (2.4) and (2.5) with Tj = iTi^E + ''"o (we will drop the e in remainder of this 
section). Since we are interested in behavior the correlation functions at late times, we will take 
Tq — )• at the end of all our calculations. 

2.3 Corrections to the Power Spectrum 

With the coupling D ;^^^^^7rcC', we expect a correction to the power spectrum Vc_ ~ A^/A;^ 
at order (jy)^~^^. The action does not depend explicitly on time and therefore we expect this 
correction to be scale invariant on general grounds. In this subsection, we will confirm this 
intuition with an explicit calculation. 

The correction to the power spectrum arises from 

bV^ = |iJ|/-2^(Ck(ro)C-k(ro) J dndT2 a(ro + in)MTo + ira)^ (2.14) 
X 7rk(ro + iri)7T_k(To + «T2)) (0-k(To + iri)C'_k(T"o + m))' , 
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where the prime on {OO)' indicates that we drop the {2tt)^ times a delta function in momentum. 
We will use the identity 

1 _ (2^)2 r(2- A) f d^k „,fc.,„,,A-2 



(x2)A - 4A-1 r(A) J (27r)4'' ^ ' ^^-^^^ 
where x is a 4-vector. The two-point function of O can be written in momentum space as 

(Ok(ro + m)0_k(ro + ir2))'= (2.16) 
(27r)2r(2-A)__,_ , ^-A.,. , ^-A 



4A-1 r(A) 



Plugging back in (2.14) the two integrals in ri and T2 can be done analytically, and we are left 
with the integral in lo to be done numerically. The correction to the power spectrum is then given 

by 



/i \ 4-2A 



SVc ^ t{A) (2.17) 

where 

^ ' r(A) y (ld2 + i)2+a ^ ^ 

(fo + i)(cD - i)e^(2fo+-A)(i _ i^)Ar(^ _ 1^ _^^(^ _ ^)) 

+(fo - 0(^1; + + icD)'^r(A - 1, -fo{Cj + i))] X 
[(fo - i)((I; - i)(l - iCj)^T{A - l,fo(w - 0) 

+(fo + + i)e*(2fo+-A)(^ ^ icD)^r(A - 1, fo{u + i))' 

where r[s,x] = t^^^e^^dt is the upper incomplete gamma function, and where we take fo = 
To/k — ?■ at the end of the calculation. In the limit fo — )• 0, t(A) is independent of k and 
therefore the power spectrum remains scale invariant. The function t(A) is plotted in figure [2j 
Notice that its numerical value is quite large, though the correction to the power spectrum is 
safely much smaller than one for a large range of values for /x. The large value is partly related to 



our conventions in equation 2.1 , which differ from free field conventions by a factor of 4tt'^ ~ 39.5. 
Further, notice the divergence as A — )• 2. This is due to the necessity of a divergent counter-term 
for the two point function for A > 2. In the language of the EFT of Inflation, the unitary gauge 
operator that provides the counter-term is (Sg^^)'^ ^M- Upon re-insertion of vr, this terms contains 
indeed the quadratic term tt^. We see that a speed of sound different from unity is generated. 
We discuss about radiative corrections and renormalization in larger detail in Appendix [Xj We 
conclude that, upon renormalization, the contribution is small even for A > 2. 

2.4 Bispectrum from (OO) 

Let us begin to explore signatures in the bispectrum. The simplest case to compute is when the 
bispectrum in n (and so in C) is induce by the power spectrum (two point function) of Cs. This 
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Figure 2: Numerically computed t{A). 



arises from a combination of two types of vertices from the following interaction Hamiltonian 
density 



n 



int 



(2.19) 



Here our operator O couples to vr both linearly and quadratically. A bispectrum can therefore 
be induced by the power spectrum of Cs. Notice that this possibility arises already in the case 
in which only the first operator in (2.19) is present. However, as we explain in Appendix |A| this 
combination cannot give rise to a large Jnl- We therefore concentrate in the combination of the 
two operators. For simplicity we work just with the operators tt and tt^, and neglect the (5j7r)^. 
The induced bispectrum takes the following form 



(2.20) 



(Cki(To)Ck2(^o)Ck3('^o) j dTidT2 a{To + iTi)^a{To + iT2f X 

'i"-ki(TO + «n)7r_k2('^o + «T-2)vi"-k3('^o + iT2))' X 
(Cki (to + «n)C-ki (to + ^^2))' + permutations . 

A detailed understanding of the shape of this bispectrum will be the focus of remainder of this 
subsection. For our analysis it will be useful to work with the equivalent expression 



B 



2-A 



h) [^H ) 8 kik2k3 
j dT2 {-iT2 - ro)^ (1 + iA;2ro^)e-^'^l^2l (1 + ik3ro^^)e-^^\-^ 
{2n) r(2 - A) 

4A-1 r(A) 



dr, (-zn - ro)^-2(l + iA;iro^)e-'=il-il (2.21) 

I 



dw e*^^i2 ^ ^2^A-2 ^ permutations 



Here we have used (2.15) to express the {(DO)' correlation function in Fourier space (and have 
also introduced more compact notation B = B{ki,k2,k3)). 
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2.4.1 The Squeezed Limit 

Let us start by analyzing the squeezed limit. Understanding the scaling behavior of the bis- 
pectrum in the squeezed limit is instructive, both because it is an important signature of these 
models and because it is possible to perform analytically. 



Before beginning to analyse the formula (2.21), let us give a bit of intuition on how a contri- 



bution in the squeezed that is larger than in single clock inflation is generated. Let us concentrate 
in the squeezed limit fci <C ^2 — ^3 • The reason why in single clock inflation there is a vanishingly 
small squeezed limit is that a long mode is locally unobservable, and so it cannot physically affect 
the correlation of two short modes. This determines the squeezed limit behavior. If we want to 
have a different squeezed limit, we need therefore to have locally observable long wavelength 
fluctuations. These are given by the long wavelength correlation of the Cs|^ This means that 
the leading contribution in the squeezed limit arises when the correlation function of the O's is 



evaluated at momentum ki (i.e. we can drop the permutations in equation 2.20). Let us therefore 
concentrate of these terms. The to integral is clearly peaked at a; 1/ti2. Furthermore, because 
the integral has support only for |ri^2| ^ (^1,2)^^ respectively, > l/^i- We therefore can 
schematically write 

B ~ j dn {-in - ro)^-2(l + iA:iro^)e-'=il^il (2.22) 

J dT2 i-iT2 - ro)^ (1 + iA:2ro^)e-'=^l-^l (1 + zfc3ro^)e-^3k.| ^-1 ^2(2-A) 

Let us consider the case A < 2 first. First consider the T2 integral. The integrand grows 
at least as fast as |t2|'^, which means that it is dominated by the largest possible values of T2: 
\t2\ ~ l/k2. At this point, the integral in ti goes as ItiI^""*^ for |ri| < 1/A;2, and as |tiP~^ for 
kil ^ ^21- The Ti is dominated by the largest possible value of ri as well: ti ~ 1/ki. Putting 
these scaling together, we obtain the squeezed limit of the three point function to take the form 

Bikuk2,k3) oc ^(^)^ forA<2. (2.23) 

As advertised, this scaling depends directly on the conformal dimension of O. 

For the case A > 2, the T2 integral behaves in the same way, being peaked at \t2\ ~ l/k2. 
Instead, the ri integral becomes a decreasing function of |ri| for |ri| > \t2\. This means that even 
the Ti integral is now peaked at |ti| ~ 1/A;2. This gives the following squeezed limit 

B{ki,k2M) oc T^f^)' forA>2. (2.24) 



k\k^ \k2 

Notice that for A > 2 the integral has a UV divergence for T12 — )• 0. We discuss this more in 
detail in Appendix |a| but here we just notice that the counterterm is tt^, which has the same 



squeezed limit as (2.24) [15]. 



*In the intuitive language described in [T^ and in [2], one can interpret the induced three-point function in the 
foUowing way. The long-wavelength vacuum fluctuations of Cs generate a long n mode by the coupling ttO, and 
affects the evolution of a short tt mode by its coupling tt^C In the example of section 



2.5 



it will instead be the 

correlated vacuum fluctuations of three Cs to generate directly a correlation between three tts though the mixing. 
The scaling in the squeezed limit can then be intuitively understood from the time-dependence of O [6]. 
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The distinctive squeezed limit that we find here and also in the next section has interesting 
consequences from the observational point of view. As noticed for the first time in [20] in the case 
non-Gaussianities of the local kind, when the squeezed limit of the bispectrum goes as A;j~^/c^^, 
the bias of dark matter halos receives a contribution that scales as relative to the standard 
bias (see [21] for a generalization to the full general relativistic setting). This is opened up 
the possibility to measure non-Gaussianities from the power spectrum of large scale structures. 
Analysis on the power spectra from current data in |22j have produced constraints comparable 
to the ones from CMB, while analysis of the bispectrum including the scale-dependent bias are 
expected to improve these limits even by about an order of magnitude [23] . Models of quasi single 
field inflation have a more general squeezed limit, so that the scale dependence of the bias goes 
as 1/A;" with 1/2 < a < 2. Non-Guassianities from our conformally coupled sector are able to 
generate a scale dependent bias in the different interval 0<a = 2 — A<1: in a sense, they allow 
us to fill the whole range. While a detection of the non-Gaussianity induced by these operators 
can be clearly within reach, the actual detection of these smaller values of a in the data seems to 
be quite hard; see [211 HSj for first forecasts using the scale dependent bias of the power spectrum 
of galaxies. It would be interesting to see if the prospects of detection will improve upon inclusion 
in the forecasts of the bispectrum of galaxies, as it was done for local non-Gaussianities in [23j . 
as well as of techniques that tend to reduce the cosmic variance [26] . 



2.4.2 The Shape and Amplitude of the Bispectrum 

In the previous subsection, we set up the calculation of the bispectrum and analyzed its squeezed 
limit analytically. Here we will compute its shape numerically for various values of A and discuss 
the observational implications. 

The bispectra that arise here are exactly scale invariant to the level of approximation we are 
considering. As a result, the shape function defined in [27j determimes the signal-to-noise of the 
bispectrum. It is given by S{x\^X2) = X2, l)a;fx|, where xi = ki/k^ and X2 = k2/k3; the 

squeezed limit corresponds to xi — ?■ and X2 — )• 1. In the examples we have discussed here, the 
shape functions scales as 

lim 5(xi,X2) oc xf-^ . (2.25) 

rEl— >-0,3:2— >1 

As a result, there is very little signal-to- noise in the squeezed limit for CFTs satisfying the 
unitarity bound A > 1, apart from the introduction of a scale dependent bias. For this reason, one 
should suspect that these models give equilateral or orthogonal type non-gaussianityj^ Evaluating 



the integrals in our expression (2.21) for bispectrum is not straightforward to do analytically 



in general, so we will compute it numerically. The two time integrals in (2.21) can be done 



analytically, and one is left to perform only the ui integral numerically A plot of the shape for 
A = 1 and 3/2 is shown in figure [3j It is difficult to discern by eye the distinction between the 
squeezed limit in this figure and the one of single clock inflation. This is an artefact of the plot 
related to the fact that the shapes starts to have the correct asymptotic squeezed limit only for 

®In section [3] we will exhibit unitary theories with more general scaling, including ones which simply generate 
local /jvL to good approximation, with support in the squeezed limit. 

^°We do not give here the result after the two time integrations, as it is just better to let your favorite 
Mathematica-like code do them for you. 
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x?j < 10~^'~^. Such a delayed onset of the asymptotic squeezed hmit makes it hard to read it in 
the plot, but we have verified it numerically. 

Next, we will employ the optimal method for comparing shapes developed in ^27|, to which 
we refer the reader for details. The cosines, defined in equation (19) of [l^, between our shape 
for various values of A and the equilateral [28], orthogonal |15] . and local templates are shown 
in Table [T} We find that the results are consistent with equilateral or orthogonal shapes for the 
range of A allowed by the unitarity bound. 

The appearance of the orthogonal shape for A ~ | is a special feature of the bispectrum 
generated by the {OO) two-point function, something which will not have a parallel in the next 
subsection when we consider the bispectrum generated by {OOO). While this paper was in 
preparation, we received the latest results from WMAP indicating a 2.45cr hint of orthogonal 
non-Gaussianity [29]. Although this is not statistically significant with current data, it is worth 
noting the appearance of the orthogonal shape (with negligible support on the other templates) 
for a particular value A ~ | of the operator dimension. 



Table 1: Cosine of shape with standard templates for operators of various dimensions. 



Dimension 


COs(S'a, 'S'equilateral) 


cos (5a, S'orthogonal) 


COs(S'A,5'local) 


A = 1 


0.95 


-0.07 


0.45 


A = f 


0.05 


0.78 


0.59 


A = i 


0.93 


-0.15 


0.45 



The value of /nl can be defined in the standard way: 

5 {QkCkCkY 
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A3 



(2.26) 



From the interaction terms leading to our general expression (2.21 ), it is easy to read off that f^i 
scales as {^/Hf^^ {fi/Hy^. We plot the value of /nlA^/ [ij^Hf'^ {ji/Hy^^ function 
of A in Fig. |4| 
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We see the expected divergence as A — )• 2, discussed above in the in the previous subsection 
and appendix A. Note that the vanishing value of /nl around A ~ 5/4 does not indicate that the 




Figure 4: /jvl 3s a function of A. 



non-gaussianity vanishes, just that it goes to zero in the equilateral limit. We see this explicitly 
from full shape for A ~ 5/4, which is similar to the orthogonal template. 

A=5/4 



5(1,1) 




Figure 5: the shape function for A = 5/4. 



Finally, let us briefly comment on the size of /nl in a simple example like (1.5). If both fx 



4 



and fl were generated at a common scale, M^,, it is natural to expect fjr = (p/M^ and fl 

(t)^/M^+^. However, for A < 2, the constraint /i < implies that {^/Hf^^ {fl/Hy^ < A^~^ 
and therefore the bispectrum generated by (OO) satisfies /nl ^ 1- Achieving large /nl is still 
possible but requires multiple scales (e.g. and M^:) or a different UV completions altogether. 
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2.5 Bispectrum from {OOO) 

In this subsection, we will consider the bispectrum of tt that is generated by TiirA = 

and the three point function of O. The leading contribution to the bispectrum is given by 



TT-ki (to + iri)7r_k2 (tq + iT2)'k-\,^ (tq + ir^))' (Oki (tq + ^rl)Ok2 (^o + iT2)Ok^ (tq + ira))' , 

where, again, B = B{ki, k2, k^) = (CkiCk2Ck3)'- As in the previous subsection, we will start by 
considering the squeezed limit before discussing the full shape. 

2.5.1 The Squeezed Limit 

The calculation of the bispectrum simplifies in the squeezed limit, ki <C k2,ks, if this limit also 
corresponds in the OPE limit of the 3 point function of O in position space: |x2— X3I ^ |xi— X2I ~ 
|xi — X3I and similarly |r2 — r3| <C |ti — r2| ~ |ri — rsl- To estimate when these limits coincide, 
let us write 



The integrals over 2:13 and X23 will receive most of their support from X23 ^ A;^ and X13 < k^ 
due to the oscillations at larger values. Similarly, the integrals over are exponentially suppressed 
unless \Ti\ < k^^. In the squeezed limit, /ci <C ^2,^3. As a result, the integral receives support 
from Ti ^ T2,T3 and xi ^ X2,xs. If these configurations dominate the integral, then we can use 
the OPE limit of the three point function to approximate the squeezed limit. 

We want to determine if the integral receives its dominant contribution in the OPE limit. To 
do so, we may assume that 2:2, 2:3, T2, Ts ~ /c^^ ~ k^^ <C k^^ as larger values arc suppressed. 
Now we will integrate over X13 at fixed ti, assuming X12 ~ X13. This assumption is reliable is the 
integral is dominated by xi ':$> X2,X3. For X13 <^ T13 the integral scales as ^^3 and therefore it is 
dominated by the largest values of X13. For a;i3 S> T13 the integral scales as xf^^^. For A < |, 
the largest values of X13 still dominate up to the cutoff, yielding X13 ~ k^^. Having performed 
the a:i3 integral, we now perform the integral over ri to find ti k^^ . As a result, when A < | 
the integrals are dominated where the OPE limit of the CFT is applicable. 

For A > |, the integral over X13 is peaked at X13 ~ ri3. Now we perform the integral ri with 
a^i3 ~ Ti3, which scales as / dri For A < 2, this integral is peaked at n ~ /c^ ^ 

where the integral is cutoff by the exponential. In this case, the OPE limit is again applicable. 
However, for A > 2, the integral is peaked around ri ~ A;^^ where the OPE limit does not 



(Cki (ro)Ck2(ro)Ck3(To) 



dTidT2dTz a{TQ + iTi)^a{To + ir2)^a(ro + its)^ x 



(2.27) 
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apply Schematically, these results imply 

Now let us compute the A < 2 case more carefully for the squeezed limit. We have established 
that the OPE limit of the three-point function is where the integral is dominated, so we can take 

^3 3 

/^gjki-Xi3+jk2-X23 
(1^X1^(^X2?, ^ 
(x?3 + r23)A(^2^ + ^|3)2 

f (n / ^r.Hr.-.o)-^(l + ....o^)e-»N) (2.32) 

i=l 

r(2 - A) r(2 - I) /" dujidL02 



x(^^) 22A-2r(A) 2A-2r(f)y (27r)2 " K+^i) (^^2 + ^:2) 2 



■2 



2 

We can evaluate the Tj integrals using 

/oo 
dr(-zr - ro)^-2(l + iA.ro^)e-'=l^l+*-^ (2.33) 
-oo 

• -(ifcro+roW+YA) 

= ^^y.^ ^ e^'^^^+'''^{k + iuj){k-iu:)^{i + kTQ)T[^-l,iTo{k + iu)] 

(k^ + w^)^ L 

-{k-iu}){k + iuj)^{-i + kTo)T[l^-l,-iTQ{k-iu))] . (2.34) 
The squeezed bispectrum is then given by 



where 



^ = -(f) T^yi2^ ("^ + ^1) (-2 + ^3)2 

X J(A;i,u;i,ro)J'(A;2,t^2,To)J'(A;3, -wi - W2,ro) ,(2.35) 

. .4 r(2-A) r(2-|) 
^^ = (2-) 22A-2r(A)2A-2r(A) (2.36) 

In the squeezed limit, the integral is dominated by L02 ~ k2 and a;i ~ ki so we may set ijji+ijj2 — s- a;2 
and ks ^ k2- The remaining A;i_2 is dependence determined in the tq — )• by rescaling Ui by A^j, 
such that the integrals are fc- independent. The squeezed limit is then determined by 

(^)""a?^3^/(A) (2.37) 



^^In a correlation function in a CFT involving more than two operators inserted at various points, the OPE 
between two operators inserted at two points is convergent as long as there is ball that contains those two points 
but no other one. 
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where the function /(A) is shown in Fig. [6] in the blue un-dotted hne. In the last passage, we 
have parametrized C = 5a(A — 1). This parametrization originates from the fact that consistency 
of the CFT under crossing symmetry requires the following numerically-found upper bound 
C < 5 a (A — 1), where a is bounded above by a number numerically close to one |30l I31j . This 
result is relatively intuitive, as A = 1 represents a free theory, and the bound just quoted implies 
that the C is continuous in the limit A — t- 1. 

The rapid growth of /(A) near A = 2 is not physical, but represents the increase in the error 
we are making by using only the leading term in the OPE. The breakdown in the OPE at A = 2 
appears as a logarithmic divergence in /(A) at A = 2. This divergence is removed if we compute 
the squeezed limit using the generally valid formula that we explain in the next section. We 
plot this function in Fig. |6] in the red dotted line, where we see that the limit A — ?■ 2 is smooth 
and the contribution remains finite. We expect that the difference between the analytic and the 
numerical functions for A —t- 1 is a numerical uncertainty. 



-1000 



-2OO0 



-3000 



-4000 



-5000 




Figure 6: Numerically computed /(A). 



2.5.2 The Shape and Amplitude of the Bispectrum 



Now let us consider the shape of the bispectrum as a function of A. As we did in Sec. 2.4.2[ 
we will numerically plot the shape function, S{xi,X2), and compare it to the standard templates 
using the cosines defined in equation (19) of [27J. We will also determine the value oi Jnl, finding 
it to be naturally substantial while at the same time bounded in an interesting way by limits on 
the size C of CFT three-point functions found using crossing symmetry. 

Computing the bispectrum analytically is not straightforward, but it can be computed nu- 
merically. To speed up the integration, one can perform some of the integrations analytically, 
as we discussion in Appendix [Bj Aside from the scaling in the squeezed limit, the shapes do not 
change dramatically as a function of A. A plot of the shape for A = 1 and 2 is shown in Fig. [7j 

The cosine, defined in |27] . between our shape for A = 2 and the equilateral [28], orthogo- 
nal [15], and local templates are shown in Table [2j The cosine varies little as a function of A for 



This is sometimes referred to as 'bootstrap' or as 'OPE associativity constraint' 
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10 .V| 10 ,v| 

Figure 7: Numerically computed shape function, S(xi,X2) evaluated for A = 2 (top) and A = 1 (bottom) 

all values consistent with the unitarity bound, A > 1. As a result, we find that our bispectrum 
is largely equilateral in shape. 



Table 2: Cosine of shape with standard templates for operators of various dimensions. 



Dimension 


COs(S'a, 'S'cquilatcral) 


cos ( 5*^ , S'orthogonal ) 


COs(S'A,5'local) 


A = 1 


0.94 


-0.11 


0.47 


A = 2 


0.90 


-0.25 


0.56 



The amplitude of the bispectrum, /tvl, clearly scales as C {n/ H) ~ since it is proportional 
to the CFT three-point function and three insertions of the mixing interaction. As discussed 
above, C is bounded as C < 5a (A — 1), where the coefficient a is bounded above by a number 
numerically close to one [301 E]. We plot the value of /^LA^/(/i/if)6-3A/ 

a as a function of A 

in Fig. [Sj Since a < 1, the curve can be thought as an upper bound to the value of fNL- We see 
that for A — )• 1, the fNL induced by {OOO) goes to zero. However, C can be of order one already 
as soon as A > 1.2, so large values of Jml easily fit within the parameter range of consistent 
theories. 

In particular, for C of order 1 the amplitude of the (nearly equilateral) non-Gaussianity gives 
us a sensitive probe of higher dimension operators such as suppressed by a high scale M^, 
which can be much larger than H. Similar remarks apply for all the cases we analyze in this 
paper. 



2.6 The Collapsed Limit of the Tri-Spectrum 

Another important signal of these models is the collapsed limit of the tri-spectrum. The behavior 
in this limit is often parameterized for local-type non-gaussanity in terms of tnl which is defined 
as 

^ _ 1 l;™ (CkiCk2Ck3Ck4) /2 oox 

4|ki+k.KoPc(^i)n(^3)n(|ki + k2|) ■ ^' ^ 

If tnl > (|/nl)^ (the I is a result of the conventions of /nl and tnl), then more than one field 
must contribute to C [211 [Ml [SI [Ml [Ml [M] ■ Such models also have important consequences in 
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/nl(A) ^f{^^m-'■*^''|a 




Figure 8: /jvl a function of A. 



large scale structure, where they lead to scale-dependent stochastic bias [381 ES]- These result 
generalize straightforwardly to tri-spectra that scale as any inverse power of | ki + k2 1 in the limit 
|ki + — (see e.g. [36J). Our models include such generalizations so we will define tnl to be 

r Urn (fclfc3)^' (Ck,Ck,Ck3Ck4)^ ,2 3ql 

"""^ 4|ki+k.Ko|ki+k2|2A.p^(A:i)P^(A:3)Pc(|ki + k2|) ' ' 

where Aj is a parameter that is determined from the leading contribution to this limit. In Section 



2.6.2, we will find that Aj is the lowest dimension appearing in the OPE of 0{x)0{Q). 

Whenever extra degrees of freedom contribute significantly to the squeezed limit of the bis- 
pectrum without contributing significantly to the power spectrum, a large contribution to t^l 
is expected [36]. Specifically, under such circumstances, one can expect to find tnl ^ (I/nl)^, 



^5- 

where f^L here is the one induced by (OOO). The basic intuition is that generating /nl requires 
three vrO mixing interactions, whereas tnl gets a leading contribution from four mixing interac- 
tions combined with the CFT four-point function. As a result, the tri-spectrum is enhanced by 
two inverse powers of the coupling controlling the mixing. In this section, we will confirm this 
intuition for our CFT examples. 

We will discuss two concrete possibilities for a large tri-spectrum. First, we will consider the 
tri-spectrum induced by the nonlinear coupling of O to vr^ that we discussed in Section 2.4, The 



full tri-spectrum is computable in this case, but we will focus on the collapsed limit. Second, 



we will consider the tri-spectrum induced by the linear mixing interaction used in Section 2.5 
Here the tri-spectrum is determined by {OOOO) and is not determined by conformal invariance 
alone. The collapsed limit can be understood in terms of the OPE of O. 

2.6.1 The (C'C')-induced Tri-Spectrum 

From the discussion in Appendix [a| it is possible that the only interaction term isH = 77^0/2/2^. 
Under these circumstances, the bispectrum would vanish, but a large tri-spectrum could be 
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generated as follows. We will define y = (CkiCk2Ck3Ck4)') and we get 

y = /i-^^MpMi7|2 / ^ / (2.40) 

X (Cki(ro)7r'_ki(n))'(Ck2('ro)7r^k2(n))'(Ck3(T-o)vr^k3(^2))'(Ck4(To)vr^k4(^2))' 



X 



(Oki+k2(Ti)0-k,-k2(r2))' 



where again {■■■)' indicates dropping the momentum conserving delta function. It is straightfor- 
ward to plug in the two point functions of tt and O to find 



I dT2 {-iT2 - ro)^ (1 + ihTo^{) (1 + iA;4To^)e-('=3+^^)l^^l (2.41) 

^ jS ^r(A)^^ / ^'^^ + |ki + ka 1 2)^-2 + permutations . 

From here, it is straightforward to compute the full tri-spectrum numerically but we will not 
show the result here. 

In order to gain more insight into the form of the tri-spectrum, we will consider the collapsed 
limit, |ki + kal — s- 0, analytically. Because the exponential suppression, ti < k'^^ and T2 < . 
Using |ki + kal ^ ki,k3, we see that for A < | the oj integral is dominated by w ~ |ki + ka]. 
Therefore, we can ignore the factor of 6*"^^^^ and compute all three integrals analytically. As a 
result, we find 

y ^ (^Ij A^^A(l + cos(vrA))r(3/2-A)r(l + A)^3±^|^, (2.42) 

We see that tri-spectrum scales as k^^~'^k^^~'^\]i.i + kap^^^ in the limit |ki + ka] — )• 0. For 
A > 3/2, the tri-spectrum scales as kl~^kl~^{Max[k3, ki])^^-'"^. Using A^ = A in (|2 39[), we see 



n \„-(ki+k2)\Ti\ 



/ - \ -2A 

that tnl oc ( j A^ and is potentially in the measurable range. Furthermore, because the 
bispectrum vanishes at tree level and radiative corrections can be small, as discussed in Appendix 
[At this can be the leading source of non-gaussanity. 

2.6.2 The (CO CO) -induced Tri-Spectrum 

Now let us repeat the calculation of the tri-spectrum (again defining 3^ = (CkiCk2Ck3Ck4)') using 
only the interaction Hint = cO. In this case, the tri-spectrum gets its leading contribu- 

tion from {OOOO) and can be written in general as 

y = /-^^Mp^|i/|2(n / ff3(_i^^^_^^)3 (Ck.(ro)vr^k,(ri))^)(0k,(ri)0k2(r2)0k3(r3)0k4(r4))(2.43) 



i=l ■ 



where again {■■■)' indicates dropping the momentum conserving delta function. Unlike the 2- 
and 3-point functions, the 4-point function in a CFT is not determined by symmetry in general. 
However, we can say something general in the collapsed limit, where ki + ka — t- 0. As with the 



20 



squeezed limit of the bispectrum, for some range of A, we expect this wiU correspond to the OPE 
hmit of the 4 point fmiction, where xi — )• X2 and X3 — )• X4. 
We will take the OPE limit first, using 

O(x)O(0) ~ ^ x^O, (2.44) 

i 

and then check the circumstances under which the corrections will be small. We can write the 
4-point function, in this limit as 



A, 



(?.45) 



where the sum runs over scalar operators of dimension Aj. We are ignoring operators with spin 
because the unitarity bounds imply that their dimensions will be > 3, and so, as we will see, they 
contribute sub-dominantly in the collapsed limit. To make the collapsed limit clear, we made the 
change of variables X14 — )• X12 + X24. In this limit, the tri-spectrum takes the form 



i=l 

d3xi2d^X24d^X34|CiPe^''i-''l2+i(ki+k2)-X24+ik3-X34 

/ ■ (^-^^^ 

Because fcj S> |ki + k2|, it is easy to see that |rj| ^k^^ ■ Therefore, r24 ~ A;2^4 and is unimportant 
for the collapsed limit (notice the difference from the squeezed limit). Therefore, the question 
of whether the OPE limit is applicable is determined entirely by the X24 integral. For 3:24 > r24 
this integral scales as X24 When 2Aj > 3, this integral is dominated by the smallest values, 
namely X24 ~ T24 and the OPE limit is not a good approximation. To use the OPE to compute 
the collapsed limit, we must have a scalar operator of dimension Aj < |. 

Following the same procedure as before, we can use equation (2.15) to rewrite this as 

i=l 

,_,_(27r)2r(2 + 4"-A)^2^^_^2 r(2-Ai) f dLOidcoiduj^ 

X 



Z^l «l Vo2A-A,-2r,'A _ AiN ; ^ ^ 



A, 



V22A-A,-2r(A-^)/ ' ' 22A»-2r(Ai) y {2tt) 



^gia;in2+ia;2r24+i<^3r34(^2 ^ ^2)A-^-2^^2 ^ ^2) A- %-2 ^| + k2|^ + Ujj)^'-'^ . 

We notice that the above time integrals simply give factors of J'{k,u;,TQ) that we defined in 
(2.33). Therefore, we can write 

/^n8-4A A^ V^|r'|2 /A A\ f duiiduJ2du}3 2 , 2nA-^-2/, 2 , 2nA-^-2 
X{\ki + k2\'^ + Ujl)^'-^J{ki, Ui, To) J{k2, u;2 - UJl, To) J{k3,u:3, To) J{ki, -UJ2 - UJ3,To) (2.47) 
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where we have defined 



.^ 1 /(2vr)2r(2 + ^ - A)n2 „ r(2-Ai) 



In the cohapsed hmit, we have ki ~ /c2 and ^3 ~ /C4. Furthermore, we see that the integral 

2' 



over a;2 is dominated by regions where 002 ^i, ^3 (when A j < | ) . The integral over 002 is 
straightforward and is given by 



— {\k, + k2\ +0^2) =|k,+k,|3-2A.2V^r(2-A.) • ('-'^^ 

The other integrals factorize into two copies of 



/ 



{k' + u;^)^- 2 uj, To)J{k, -co, to) = tTTTv ^0 • (2-50) 



(27r)' ' ' ^v'^,^,-u;c. ^i+A, 



Here we have dropped the /c-dependence of ^(A, Aj), as it vanishes in the tq — )• limit. Therefore, 
the tri-spectrum in the collapsed limit is given by 

3 



where ^(A, Aj) = ^^^^,^^^k{A, Ai)g{A, Ai)"^ . The function ^(A, A) is plotted in Fig. 
a given value of A, the dependence of g{A,Ai) on the value of Aj is relatively weak until we 
reach the regime Aj — t- 3/2 where the OPE approximation (and hence this calculation) breaks 
down. This origin of this breakdown is identical to the one we found in Section 2.5.1| when 



For 



calculating the squeezed limit of the bispectrum. Notice that the collapsed limit is dominated 
by the contribution of the Cs with the smallest Aj that enters in the OPE. So, whenever the 
OPE includes operators with Aj < 3/2, we obtain the leading collapsed limit. This additionally 
justifies our procedure to neglect operators with spin in the OPE. 



«(A, A) 



^ 1 ^ ^ ^ ^ 1 ^ ^ ^ ^ L ^ 

1.1 1.2 1.3 1.4 1.5 

Figure 9: Numerically computed g{A,A). 
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Finally, let us comment on the observability of the trispectrum. The signal-to-noise in the 
trispectrum scales as (C^)/(C^)^ ^-^id therefore the ratio of signal-to-noise in the trispectrum 
relative to the (C'C'O)-bispectrum is given by 



(5/iV)4 IQIVm 



{S/N)3 C \H 



2-A 



(2.52) 



For the trispectrum to give the dominant source of non-gaussianity, one requires that \Ci\^ ^> C. 
This arises naturally in CFTs which contain an approximate Z2 symmetry under which O — )• 
—O and Oi — )• Oi. When this symmetry is exact, C = and [S/N)^ vanishes (ignoring the 
(CC)-bispectrum for the moment). This is compatible with the bounds that come from OPE 
associativity [31] , which roughly bound Ci to be smaller than order one unless A or A' are within 



order 10% of 1 



3 Time Dependent Examples 

In the conformally coupled CFT of the previous section, the unitarity bound A > 1 plays an 
important role, suppressing the squeezed limit of the bispectrum in favor of approximately equilat- 
eral/orthogonal non-Gaussianity. One might naively conclude that this bound generally restricts 
the phenomenology arising from a CFT coupled to inflation. However, in the presence of time 
dependent backgrounds, infrared physics and unitarity bounds are modified [7j. Since the infla- 
tionary background is time dependent, it is important to consider this more general case. In this 



section, we will see that, because of this scaling, the squeezed limit may take the form (2.29) 
but with a larger range of powers A, including A ~ which gives local non-Gaussianity. To 
establish this it suffices to exhibit a controlled class of theories exhibiting this behavior; we will 
leave a more general analysis to future work. (Non-conformal coupling of the CFT to the de 
Sitter curvature is a natural case for further study, for example.) 

In [7| the relevant effects of time-dependent couplings were analyzed in detail for a particular 
class of calculable examples, the renormalization group flows generated by adding a product of 
two single trace scalar operators O to the Lagrangian of a large-N CFT. In the static version 
of this theory, these double trace flows connect a CFT in which an operator 0+ has scaling 
dimension A+ to one in which this operator has been replaced by an operator O- with scaling 
dimension A_ = 4 — A+. For the range of operator dimensions 1 < A+ < 2, the double-trace 
deformation by j 0\ is relevant, and the dimension of the scalar operator flows from A+ in the 
ultraviolet to A_ in the infrared. For higher A^., the deformation is irrelevant and the dimension 
of the scalar operator becomes A_ in the UV; for A+ > 3 one has A_ < 1 and this introduces 
unitarity violation at high energies. These statements are straightforward to derive in the large-N 
limit, as we will explain shortly in adapting the calculation to our application. 

As in [7], we will be concerned with the generalization to the case where the double trace 
coupling is time dependent: in particular, let us deform the CFT action by a power-law coupling 
of the form / Ai^'^O^ where t = log(ffr). For sufficiently large k, specifically k > A-|_ — 2, 



^''A very approximate fitting formula extrapolated from the Fig. 12 of [31] is Ci < (5A'/2 — 3/2)(A — 1), valid 
for A < 2 and A' < 2. 
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this renders the deformation relevant, and the theory flows to a new non-conformal theory in the 
infrared. In the leading large-N limit, one can calculate the two point function of the operator 
to which O flows in this new theory, finding a result proportional to 1/ distance^' times powers 
of the time-dependent coupling. In the next subsection, we will recover this result in our case in 
the process of generalizing the analysis to the inflationary background geometry. This will give 
a similar result for the two point function of the scalar operator, but now with an additional 
factor of (rr')^" reflecting the de Sitter redshifting. Because of that, up to logs we will recover 



the behavior (2.29) but with A replaced by A_. This realizes a larger range of powers of k^/ks, 
including the local shape for A_ = 0. 

Power-law time dependent couplings may be natural in inflation given couplings of the inflaton 
to other fields. In say rri^ip^ inflation, the inflaton rolls away from its initial value like (j) — (po oc t, 
so a linear coupling to it would naturally produce a coupling with k = 1 and so on. For fi^ip 
inflation we have oc t^/^ and one could similarly get an order- 1 value of k if powers of (p couple 
in. That said, we will not engage in serious model-building in this work; the following scenario 
is just meant to establish that strongly coupled fields can introduce a wide range of scalings in 
the squeezed limit, including the local shape. 

3.1 Modular example 

In this section, we will exhibit a concrete example which is essentially a hybrid of ^[6jj and [7]. 
The model includes a massive weakly interacting scalar field r] which mixes linearly both with 
an operator O of a large-N CFT of dimension A_|- and with the time-derivative of the scalar 
perturbation tTc- Specifically, we start from the action 

S = l^,[l{H'rHdvf-mW)+9o{H-Hog{Hr)r^O^+si,i^ (3.1) 

+ I^AHrp<V + l^v'}+S.nfi (3.2) 

with TTc the canonically normalized scalar perturbation. The time-dependent coupling g(T) = 
gQ{H~^ log^Hr))'^ = gof^ here could come from an interaction term such as J((^ — {pQ)'^r]0 
between the inflaton ip and the CFT sector. As in [6] the mixing term J pTTci] could come from 
a coupling J (dip)'^r] with one factor of d(p evaluated on the background and the other on the 
perturbation. Finally, we have included an interaction term J ^rf which will generate a simple 
contribution to the non-Gaussianity. (The CFT three-point function is down by a factor of 
so we will neglect it here but it would also be straightforward to include.) 

The top line generates a time-dependent double trace flow in the strongly-coupled sector, 
relevant for k > A_|_ — 2 as in [7j. In the regime where the mass term for rj dominates over its 
kinetic term, integrating out ry produces an operator relation 

V = (3.3) 

implying a double-trace deformation ~ J {g"^ / 2m?)0'^ . 

The two-point function for r/ (equivalently O) is given in the large-N limit, and to zeroth order 
in the cubic interaction p/ by summing up the diagrams generated by linear mixing between 
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rj and O. In flat spacetime - applicable in our case for scales well within the horizon - this gives 

ivixHx')) = ^—^ (3.4) 

where ly = A_)_ — 2 = 2 — A_ and 

It was shown in [7] how unitarity works out in the infrared in this theory. Using the technique 
in [40j, one finds that unitarity holds - it requires the positivity of the product of the —1 in the 
numerator of (3.4) and a factor (A_ — 1). The consistency of the theory at long distances is 
not a surprise; it is to be expected in a theory like this which is well-defined at shorter distance 
scales. We will now derive the corresponding result in the inflationary background of interest 
here, which will lead to the two-point function 

{t]{t, x)ry(r , x )) - ^^^_^^2(iog ^^),.(iog ^^z). [(^ _ ^,)2 + ^/)2]A- " (^'^^ 

This is obtained as follows. First, we compute the correction ~ J^, Vx{C'xOx')rix' to the effective 
action in S^ff for r] which arises through its mixing with O, giving 

f drdx I' dr'dx' r]{T,x)gl{H~^ log HtY{Ht)^+{H-^ log HT'Y{HT')^+r]{T',x') 

(3.7) 

where 5*0 is the original effective action for rj without the coupling to O. The trick is to note that 
when the double-trace term is relevant (i.e. when k > — 2), the second term in S^ff (the one 
generated by {OO)) will dominate in the infrared two-point function; in |7] this was explicitly 
demonstrated by bounding the corrections to this approximation. Therefore, to compute the two 
point function we need to just invert that term. To do that, we can first absorb the r and r' 
dependences into a new variable 

r? = g^{H-' log HtY{Ht)''+S = go{H-' log H tT {Ht)-^- (3.8) 

The second term is then simply 

drdfdr'df ^ir,x)f,{r',x') 

[{t - t'Y + {x - x'Y]^+ 

and we can invert it to obtain the two point function for fj exactly as in Poincare invariant 
Minkowski space CFT. (Specifically, one can go to momentum space where {OO) is ~ k'^^+~'^ 
and invert this to get k^~'^^+ = k'^^~~'^ times appropriate constant factors j7].) This gives the 
usual result for double trace deformations - a flow between A+ and A_ : 

(W> = ^,,,,A (3-10) 

c^c-^ [(r - T' f + {x- x'yr^ 



Finally, putting back the r dependence from (3.8), we obtain the claimed result (3.6). 
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The correction we generated to the power spectrum in this theory behaves paremetrically as 



h) 



1 



9l 



log'^ (H/k) log'' (H/k) 



(3.11) 



As we wiU see, the non-Gaussianity is determined by similar factors, but is enhanced by a factor 



of A; 



10^. We will focus on the case where AVt < Vr, for which the leading effect of the 



coupling to our strongly coupled sector is in the non-Gaussianity. 

As already mentioned, we will take the J ^rf term as the leading source of the non-Gaussianity. 
This leads to a calculation very similar to that arising in the case of additional weakly coupled 
fields [5] [6], and we will compare the two as we go. Schematically the leading contribution to the 
bispectrum is 



B{ki,k2,k-i) 




dT- 

P Griikj ; Tj , T)dr^ Gt,^C {kj ; Tj , 0) 



(3.12) 



Here Ctt^^ = GqMp\J 2\H\/ H = G(^H/A(^ is proportional to the two-point Green's function of 
the scalar perturbation C, and 



dr^G(;{kj;Tj,0) = dr^ < Ty (1 + kj\Tj\)e 



-kj\Tj\ 



(3.13) 



(working in Euclidean signature). In (3.12), Gr^ is the two-point Green's function of rj, given by 



Fourier transforming (3.6): 



25/2^3/2 ^2A_+2k k]^- 3(^^.)A_ . |^ _ ^. 



Grj{kj', T, Tj 



(3.14) 



c,c_, r(A_)<7g \k,{T - T,)\^-~yHog'^iHT)log^{HTj) 
Before continuing, let us remark on the combination 

go = Vr(A_)5o (3.15) 

which appears here since we will be particularly interested in the regime A_ — )• 0, and the F 
function has a pole there. This effect is also evident in the position-space propagator, where 
l/\x — x'p'^- = 1 — 2A_ log \x — x'\ + . . . as A_ — )• 0; the two point function of r] only depends 
on X — x' through the log piece proportional to A_ . The term independent of x — x' will not 
contribute to our three point function at finite spatial momentum k, so the leading term we need 
in the two point function of rj is the log term proportional to A_. If we take to be finite as 
A_ — )• 0, then the two-point function of ry, including the effects of its mixing with O, stays finite. 
It is this degree of freedom which couples linearly to the inflationary perturbation, and it seems 
natural to keep its propagator of order 1. This regime is consistent with the flow we are working 
with [7j occurring within the horizon. To see this, one can consult equations C.1-C.3 of [7j, whose 
leading x — x' dependent term depends on go and A_ through the combination ^o- 



From these expressions (3.13) and ( 3. 14 ) , we can now determine the leading contribution to the 
integrals in (3.12) in the squeezed limit k2 ^ k^ ^ ki. For the short modes (3.13) exponentially 
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suppresses the integral for r2,3 > 1/^2- As a result, the Bessel function i^A_ -3/2(^1''" ~ Tj\) 
suppresses the integral for r > 1/ A;2 . 



Next, consider the opposite limit of small kr in our integral. The two-point function (3.14) is 



different from the two-point function of massive fields, as is easy to see in position space where 



the latter is a hypergeometric function whereas ours is the simpler function (3.6). However, they 
behave similarly at small kr and kr', where for A_ < 3/2 ( |3.14 ) becomes 



^2A_+2Kr(3/2 - A_) {kjT)^" {kjTj)^ 



(3.16) 



Up to the logs, this is just like the corresponding expression for weakly coupled fields of mass rh 
(quasi-single-field inflation), with the identiflcation 



A_ ^ 3/2 - a/9/4 -m2/i?2 



(3.17) 



The Tj integrals in (3.12) are manifestly dominated by the largest value tj ~ l/fcj allowed by the 



exponentials. As in [5j and the example of the previous section, the small r behavior does not 
lead to large effects; again, this can be seen from the contour prescription |18j . 

Given that, the integral is well approximated by a saddle point with r ~ r2_3 

n 



kg^ and 



kr ■ Plugging this into (3.12) it is easy to read off the scaling and amplitude in the 



squeezed limit 



B{ki,k2,k3) 



with 



TNL 



fNL^l 
1.3 p 



^n3 

H 



kL\ 
ks) 



(3.18) 



50 



(3.19) 



up to logarithmic factors. Here we have dropped factors that are order 1 in the regime A_ <C 1 
which approaches the local shape. 




Figure 10: Numerically computed shape function, S'{xi,a;2) evaluated for A- = and k = 1. 



The shape can be computed numerically starting from (3.12). The case of A_ = and 



K = 1 is shown in Fig. [10} The bispectrum is not scale invariant with an effective /nl that 
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grows logarithmically at larger scales (smaller k). As a result, shape function is enhanced in 
the squeezed limit by additional factors of l/logX2, which represent the stronger interaction 
at large scales. Numerically, we find significant overlap with the local shape for all A_ < 1 
(cos(5a_ , •S'locai) ^ 0.8). As discussed around Fig. [6| the numerical error appears to grow as 
we approach the squeezed limit, which is also where the signal-to-noise is dominated. For this 
reason, the specific values of the cosine with the local shape are likely not meaningful. 

Let us compare our amplitude to the case of an interaction J tTcCt between the inflaton per- 
turbation and a weakly interacting field a of mass m. In the weakly coupled case, the local 
shape arises for massless scalars, m — )• 0. But the cubic interaction - which here produces the 
non-Gaussianity - would generate a mass m ~ /x, unless it is tuned away. Without such tuning, 
the amplitude is proportional to the mass through its dependence on /x. In our case, although 
we have the fir]"^ coupling, we do not need to take the mass m of ?? to zero; instead we get the 
local shape from A_ <C 1 rather than m <C -ff . Our amplitude Jnl can easily be large for small 
A_, with the flow we have used fitting well inside the Hubble horizon. Finally, we also note 
that the calculation of the two point function in this example requires a large-N limit to control 
as in [7], which may be regarded as a small (perhaps ten percent) tune we have introduced for 
calculational convenience. 

As a final comment on the amplitude of the non-Gaussianity, let us comment briefly on 



the sensitivity this gives to higher dimension operators in the theory (3.1). This example is 



more natural than the simple case (1.5) discussed in the introduction in that the latter involved 
a relevant operator, whereas here we may consider a nearly marginal operator C'A+Ki4 in the 
UV, with a relevant flow induced by the time dependent background. Let us make a quick 
estimate of the precision measurements available via /ml in this example. Focusing on the case 
A+ ~ 4, A_ w 0, K ~ 2 and assuming a common scale M^, suppressing the higher dimension 
operators in the UV theory, we have couplings of the form 



4^r'r,0^(34^^^ (3.20) 



M3 M, 
In terms of the parameters used above, we have 

Qo ~ ~ ~ and n ~ ~ (3.21) 

where we used the slow-roll relation A^ ~ H"^ /cj). In order for the double-trace flow to fit within 
our Hubble patch, we need = H/eq with eo < 1- This implies 



A2 



(3.22) 



Ml eo 

with 

To get a rough estimate of the sensitivity here, a null result /nl ^ ^ would probe the regime 
H and eo ~ 10"^ with M ~ 



f^,-A-^S^('^)et^Afmel (3.23) 
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4 Discussion and Future Directions 



In this work we have analyzed the non-Gaussian corrections to the inflationary perturbations 
arising from a hnear couphng to a scalar operator O in a strongly coupled theory. Specifically, 
we analyzed two simple, calculable cases: a conformally coupled CFT (at least near the Hubble 
scale), and a particular type of time-dependent deformation thereof. We found characteristic 
scaling behavior of the bispectrum near the squeezed limit, going like ~ {ki/kgy times the local 
shape. This behavior is similar to that generated by weakly coupled massive fields [5] [6], but with 
the exponent 7 in our case depending on the dimension of O. Our two examples together exhibit 
a large range of exponents 7, and a range of shapes, including equilateral and orthogonal for large 
7 and the local shape for 7 — )• 0. In the nearly-equilateral case, it is at least partially degenerate 
observationally with single field infiation (e.g. \i2\ [T3l [T^ [T5]). As in that case, the amplitude 
can naturally be large and gives us an observational probe of higher dimension operators. In 
particular, we have seen here that non-Gaussianity provides a precision test for additional sectors 
of fields coupled through higher-dimension operators suppressed by a scale M^, ^> and either 
a detection or a null result would be very informative. 

There are many interesting generalizations. Here we considered perturbative couplings be- 
tween the inflationary and strongy coupled field theory sectors, leading to small corrections to the 
power spectrum along with relatively substantial contributions to the non-Gaussianity. It would 
be interesting to consider stronger mixing interactions, requiring resummation to determine the 
power spectrum (something which simplifies somewhat if one considers a large N limit). Here 
also we focused on scalar operators, but higher spin operators could lead to their own distinctive 
effects, generalizing |41| to the case of strongly coupled fields. Another interesting direction is to 
understand more systematically how non-conformal couplings affect the bispectrum in interact- 
ing theories, a question amenable to perturbative quantum field theory calculations of anomalous 
dimensions in some interesting limits. 

In another direction, UV complete mechanisms for inflation studied thus far often involve 
strongly coupled sectors which play a key role in producing dynamically the small scales required 
in the inflaton effective action. It will be interesting to determine the implied couplings between 
their operators O and the inflaton perturbations in such examples, a potential source of new 
signatures or constraints. Related to this, it would be worthwhile to assess more systematically 
the level of Wilsonian naturalness in various models of multifleld perturbations, in preparation 
for the observational results which will determine its actual natural-ness. 
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A Radiative Corrections 

In this appendix we discuss radiative corrections to the effective action of the models considered 



in Section [2| As discussed in the main text, the relative sizes of the couplings in (A.l) determine 
which correlation functions of our operators O generate the dominant contribution to the bis- 
pectrum. Because the size of the couplings is also important for radiative corrections, we must 
evaluate the radiative stability for the different cases separately, and our analysis in this appendix 
will follow the order we take in the main text. In all cases, we will ask what is the natural value 
of /nl from the tree level Hamiltonian and compare to those generated by loops. 

A.l Case 1: Bispectrum determined by (OO) 

As in Section 12.41 we will consider the interaction Hamiltonian 



^int = |/i^-^Mpi|i?|i/2(2^ - a^^a^7r)0 + iMp2 |ij|/i-^(-2^ + df.nd^'TTfO . (A.l) 

Contribution of the second operator gives the dominant contribution to the bispectrum when 
^ fJ-^-^iMpiH)^/'^, which is the regime of interest in this subsection. Because the second 
operator is always irrelevant, the interaction becomes strongly coupled (i.e. perturbative unitarity 
breaks down) at the scale 

A[/ ~ /i . (A.2) 

Another scale in the problem is ~ Mp\H\^/'^ , the time translation symmetry breaking scale. 
Let us further assume (O^) is negligible, so that the dominant contribution to /nl is from 

(C') ~ (C^c)V'-V^(O') (A.3) 

which gives 

Clearly this can be large. Let us first consider the renormalization of the term in the Lagrangian 
proportional to O. This is generated by a loop of vr's to give 

m^^- D . (A.5) 

This is a relevant operator only for A < 4. In this case, we have to impose that 



a4 / r, \ 



^S/f"-^ ^ A<h{^] . (A.6) 



This can be clearly satisfied for fl ^ H while at the same time having A ^ H. However, if we 
push A to be as high as the Ajj, we find that the condition becomes fl < H, which cannot be 
accepted. This means that in order for having large effects from this second operator, we must 
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have A <C A.u ■ This requires an additional scale beyond A^ to arise in the UV completion of our 
model. One possibility is that that role is played by ~ Mpi|iif|^/^, which is also a physical 
scale in the system. 

Another term which is generated by renormalization is 

, /ux2-A(M2|ij|)3/2 

^ (x) A2-A^A (-2^ + ^M^^"^)' ' 

which generates a contribution to /nl of order 

f(2) ^/X\2-A/A\AMp^^Hf^ (,w//^4-2A 
~' .a) Kfi) A2 ^^^^\A. 



c(2) ^ .(1) 



For A > 2 we have f^v^ ^ /nlj which tells us that in the case A > 2 non-gaussianities from the 



second operator in (A.l) are subdominant to those induced by radiative corrections, so they do 
not dominate in technically natural theories. 

This situation will not extend to the case in the next subsection. There we find theories with 
A < 8/3 for which the theory is technically natural and the leading signal is from the three-point 
function of the conformal sector, provided the latter has sufficiently large amplitude C. For 
C ~ 0, we will find technically natural theories with large non-Gaussianities only for A < 2, and 
these will require a crossover to a UV completion at a scale below Ajj. 

A. 2 Case 2: Bispectrum determined by (OOO) 



As in Section |2.5| when {OOO) gives the dominant contribution to the bispectrum, we only 
require the interaction 

■Hint = /i^-^Mpi|i?|i/2(_2^ + d^,TTd''TT)0 . (A.9) 

For A < 2, the leading operator is relevant, and therefore, for control we require fi < H. On 
the other hand, when A > 2, the leading term is irrelevant and therefore > H is required for 
control. 



A. 2.1 Relevant deformation: A < 2 



When A < 2 only the second operator in (A.9) is irrelevant. The strong coupling scale associated 
to this interaction is therefore set by 



(Mp^il^l)^ 

2-A 



Because n < H,we see that A^ > ~ Mpi|ij|^/2. 

Schematically, the three-point function receives two tree level contributions: 



,,6-3A 



4-2A 



Mpi|F|V2 



{O' 



(A.IO) 



(A.ll) 
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Assuming the absence of UV divergences when we take these expectation values for A < 2 (a 
fact we will recover shortly) the the natural values of /nl are given by 

C ~ (A.12) 

where C is the coefficient of (O^). Here we see that only has the possibility of being large. 

Now let us look at which operators are generated by loops or by CFT dynamics, given a hard 
cutoff at the scale A. We must have H <^ A < Ajj. Schematically, we expect to generate the 
following operators 

.4 „4-2Aj^2 

^ ^'"^ Mi|g|V2 ^ + ^ A4-A {-'^^ + d.^d'^fO (A.14) 

The most relevant term is the first one and is dangerous because it can break the conformal 
symmetry. In order for this term to be negligible, the effective scale of the breaking must be 
lower than Hubble. This can be achieved, provided that 

<^H^-^ ^ A«h(^Y^\ (A.15) 



This can be satisfied for H <^ A <^ Ajj- However, we cannot push A all the way up to Au 
because in this case the constraint would become Au <^ H, which is not acceptable. We have 
not explicitly UV completed our theory; it may happen in such a completion that additional 
operators appear suppressed by the lower scale /,r ~ Mp\\H\^^'^, or at some other scale A that 
depends on the dynamics of additional sectors. We will not engage in serious model building 



here but note that our low energy model as written only pertains below the scale A <C Ajj. 

By looking at the quadratic term tt^, we see that a speed of sound slightly different from one 
is generated 

Now let us look at the contributions to /nl from the second through fourth terms in ( A.14| ) : 

& ~ c(^)'-'"(f)'-\-«/™. (A.17) 



We see that gives the dominant contribution to the bispectrum. This means that for A < 2 
the theory is technically natural. 



^something to potentially return to depending on the outcome of measurements of /n 
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A. 2. 2 Irrelevant Deformation : A > 2 



To distinguish the two cases, we will redefine ^ — t- /. Let us identify the strong coupling scale. 
Since there are two irrelevant operators, the unitary bound is given by the smallest scale sup- 
pressing these two operators. We have 

ku = Min[/,/i-2/A(Mp2if)V(2A)] ^ Min[Ac7,i,A^,2] • (A.20) 

Notice that since we must have Ku S> this implies f ^ H. Assuming no UV divergences, 
there are again two contributions to the bispectrum with 

/ill - C{j) (A.21) 



The similarity to the first case is expected here because nothing has changed, at this level. The 
differences will be more apparent when we discuss radiative corrections. 

We will again look at renormalization of various operators assuming a hard cutoff A, with 
if <C A < Ajj. As expected, we will renormalize the same operators as the previous section: 

^int'- = fL Mpifff + g /25lKliA (-2^ + d.-d'^-fO (A.23) 



A2A-4 (M2|ifn3/2 



For A > 4, the first term is irrelevant to does not introduce a constraint. However, for A < 4, 
we require 

In the case in which the unitarity bound is given by Aij^2 (that is for / > (HM^^y^^ ~ HA^ ^^^), 
A cannot be pushed all the way to A[/^2 5 as this constraint would require Au^2 ^ H. However, 

— 1/2 

when the unitarity bound is A[/^i (that is when / < HA^ ), it is possible to push A to Au^i by 
imposing the stronger constraint H <^ f <^ HA^ ^^^^ ^\ This feature of the model, that is the 
fact that we do not induce the most relevant operator even in the case where the cutoff is as high 
as the unitarity bound, can be understood in terms of an approximate vr — )• — vr, O — )• —O Z2 
symmetry that is softly broken by the operator {d7r)'^0. A very similar Z2 approximate symmetry 
is what allows to have a leading trispectrum in single field inflation [42j . 

The correction to the power spectrum scales as (A//)^^~^, which is small if A < /, which 
is always the case if / ^ (HM^^)^^^. A speed of sound different from one is generated. We 
also require A > H, which is always consistent with this constraint. Again, let us look at the 
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contributions to /nl from the second to fourth term in (A. 23) : 



(A.26) 



1 , 2, A2^-4 

g2V-^ ^s) pA-4 

(5) Mpi^A3A-6 Mpi|ff|V^ g3A-6 ^ 3A-8 
^NL ~ ^2 73A^ i72 J6A^[h) ' ^^■^^) 

For A > I we have fi^l > f^. For A > 4, f^^ > /jjj^. This means that in the case A > 8/3 
the signature from the three-point function of conformal operators is not the leading one for 
technicaUy natural theories. 

A. 3 Case 3: Trispectrum determined by (OO) 



If we start with the Lagrangian (A.l) with /i = 0, we have that the leading signal in the trispec- 



trum induced by (OO) with a size tnl^^ ~ if^/H)"^'^ . Naively, the operator has a Z2 symmetry 
under which vr — t- — vr that forbids the appearance of any odd power of vr. This naively forbids 
the presence of operators that would induce a non-vanishing bisectrum. However, the non-linear 
realization of time-diffs. imposes that the operator {6g^^)'^0 D vr^O also contains cubic terms 
in vr. Under a loop, these terms will generate terms of the form nO. If A < Ajj ~ /i is the cutoff 
of the loops, we have 

^Lt- D ^.^,~,^{HM^,,y/\0 , ~ ^ . (A.29) 

Let us see under which conditions on the cutoff the induced bispectrum is smaller than the 
trispectrum. We have 

AOOO) C ^/irad.^3(2-A) C A^^ f 
J] " 



Imposing the signal from both to be smaller than the one on the trispectrum implies respectively: 

We see that this leaves large room for A ^ H with the trispectrum being naturally the leading 
signal. Notice that already imposing the signal from the trispectrum to be detectable, tnlAi^ > 1, 
implies /i^ ~ Afj < H'^/A^^^. So the limits above represent a mildly stronger constraint. 

B Details of the Shape Calculation 



In section Section 2.5.2 we discussed the shape of the bispectrum computed numerically from (2.28) 
To simplify the numerical calculation, one can perform several integrals analytically before per- 
forming the numerical integration. This is important because naive numerical integration leads 
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to a 6-dimensional integral that is hard to evaluate. In this appendix we will explain which 
integrations were performed analytically. 

First we need to compute the Fourier transform of the three point function at different times 
and different fc's. Let us work in Minkowski space first. We know that to write it in de Sitter 
space we simply need to multiply by the relevant conformal factors at the end. Let us start by 
going to 4d Fourier space: 



(Ok.(ri)Ok.(r2)Ok3(r3)) = / e^^--- {0,.0,.0,.) , (B.l) 

where all correlation functions are meant to be anti-time-ordered. 

We know how to write the 3-point function in real space. So, let us do a 4d Fourier transform: 

f duJidU2dW3 /■ 4 4 4 ,j^fcM,^ 1 1 1 

/ 27r 27r 27r / X3 e o . {d.zj 

J J ^2 ^2 ^2 



We know how to do the Fourier transform of the 2-point function see eq (2.15), so let 

us write it in this way. It looks like we are adding a lot of integrals, but, as we will see, many of 
them can be done analytically. We obtain 

f duJi duJ2 du)-3, -i-s- u^ r- I ,4 ,4 ,4 -iTfe^'x ( (^Tt)^ r(2 - |^ ) \ 

= — — — e y^Wx^d^xae ^^M.. \^-^^^^) - 

f d^ d^ ,[p.,^^^+pM,^3^+p.,3^^] /2N A/2-2 /2^ A/2-2 /2^ A/2-2 

J (27r)4 (27r)4 (27r)4 ^^^^ ^ 

The Xj-integrals lead to three 5-four functions of the form 

(27r)^5(^)(^ fcf ) (2^)^5(^)(p^ - (pi; + fcf )) (2vr)^5W(p^ - (p^, - kt;)) . (B.4) 

Notice that 

(2^)V^)(5;^f) = (2^)¥3)(j;k,) {27t)5{Y,u^,) . (B.5) 

i i 

The first term is the usual 3-delta function of spatial-momentum conservation. We can drop it 
by adding a ' to our correlation function. Now the integral in pi and p2 can be done saturating 
the (5-function, and we are left with 

3 



.42-1 r(f ) 

d^P3 (, , , .2\^/2-2,, .2AA/2-2/ 2aA/2-2 

{iP3 + ki) ) {{P3 - ks) ) [Ps) 

The integral in d'^p^ can be done after inserting two Feynman parameters. Let us use that (see 
for example App. F of [43J): 

A-'^B-^C-^ = ^^,^" + ^+7^, dn dv X (B.7) 
n"-i(l - uf-^v^'^^'^il - vy-^ [uvA + (1 - u)vB + (1 - i;)C]-°-^-^ . 
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We then get: 

= c /^M^^ ^^^^^^ ^,.,.7(^ 

7 27r 2^ 27r ^ ' " r(A/2) ^ r(2 - f )3 ^ 







where 



M2(u,t;,{A;j}) = (nt!(l - v)kl + uv'^{l - u)kj + v{l - u){l - v)kl) , kf = k'^h^^ = Jf + kf, 

(B.9) 

and where we have shifted the variable of integration to complete the square and used that 



^ i.2 _ r,2 _ p 



The P2, integrals is equal to 



y (27r)4 ^ ; (47r)2r(6 - ^) ^ ^ 



We obtain: 



(Ok,(ri)Ok,(r2)Ok3(r3)> = ^ j ^^^)^ ^ r(A/2) 



(47r)2r(2- A/2)3 

The V integral can be done analytically (by Mathematica, not by us!), finally obtaining: 

(Ok.(ri)Ok,(r2)Ok3(r3))' = -C- ^ \l ^' esc (^A/2) T (4 - 3A/2) T (2 - A/2f 



((1 - u)n)^~t (n [kj + ujI- Loi) + ^^(-(u - 1)) + wf) ^ ^ x 

F (a_^ ^.g-l (n - l)n (fe2 + a;2) \ 

' M 2 '2'^'' {u-l)kl + {u-l)u;l-u{kf+iof) r ^ ^ 



Now we are ready to deal with our expression for the bispectrum (2.27). We notice that the 
''"1,2,3 integrals can be done analytically. Then we can do analytically also the integral over one of 
the w's, as it is just given by the saturation of the delta function. We do not give the result here 
directly, as it is not illuminating. It is just better to ask Mathematica to do it. At this point we 
are left with three integrals to do: one over u, and two over the remaining two w's. By rescaling 
the variables, it is easy to see that the bispectrum scales as k~^, signalling its scale invariance. 
We can do the resulting integrals only numerically, and Mathematica does them in a few seconds. 
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